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Abstract 

A serious problem with the Schwinger-Dyson approach to dynamical mass gener- 
ation in QED3 at finite temperature is that the contribution from the transverse 
part of the photon propagator, in the landau gauge, leads to infrared divergences 
in both the mass function T, and the wavefunction renormalization A. We show 
how, by using a simple choice of vertex ansatz and a choice of non-local gauge 
(the 'D-gauge'), both A and S can be made i-r finite. We formulate an equation 
for the physical mass A4 = '^/A, and show that it reduces to the corresponding 
equation obtained in the Landau gauge in the constant physical mass approxima- 
tion j\4 = A4{0,7rT) (which is finite). Therefore at finite temperature, we are 
able to justify a 'constant' mass approximation for A4, and show that the value 
q£ J, _ ^-^P^'^) remains close to the value obtained in previous calculations which 
included retardation. 
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1. Introduction 



A considerable amount of work has been done on investigating fermion 
mass generation in QED3 both at zero, [l]-[9], and finite temperature [10]- 
[16]. 

At finite temperature, one finds a critical temperature, T^, above which 
the generated fermion mass E is always zero. An important parameter in 
the theory is r, which is the ratio of 2E (evaluated at zero 3-momentum 
and temperature) to Tc. This quantity has been calculated in several finite 
temperature calculations [10]- [16], and turns out to have a value significantly 
higher than that found in 4-fermion models of B.C.S type, where r is ap- 
proximately equal to 3.5. Thus r can serve as an indication of the pairing 
mechanism (gauge interaction or 4-fermion) in possible applications of QED^ 
to planar superconductivity. 

So far, evidence has been found to support the idea that retardation 
causes a significant reduction in r, unlike other refinements which have been 
attempted. In [10], where only the longitudinal part of the non-retarded 
propagator was used in the constant mass approximation (i.c assuming the 
mass is independent of 3-momentum), r was found to be roughly around 
10. When [10] was refined by the introduction of momentum dependence 
in the mass [11], there was seen to be little change in r. In [12] it was 
found again that the value of r stayed roughly the same when wave function 
renormalization was introduced. However, when retardation was introduced, 
in the calculations of [13] and [14], this value fell to a value of r ~ 6. On the 
introduction of frequency dependence in addition to retardation, [14] , it was 
seen that r only changed by a small amount. Both [13] and [14], however, 
suffered from (different) major technical problems. 

In [13] (the first calculation to be done taking account of retardation), 
where the real time formalism was used, the main difficulty was the non- 
analyticity of the vacuum polarization at the origin in /c-space. In [13] such 
behaviour was only able to be treated crudely. 

In [14] the Matsubara formalism was used. The major problem for such 
calculations in the Landau gauge is that there is an i-r divergence in the 
equation for S which comes from the transverse part of the photon propaga- 
tor; so [14] neglected this term. In [16] it was shown that the introduction of 
a fermion mass in the calculation of IIi,^, does not get rid of this divergence. 

In the real time formalism a more refined calculation has recently been 
attempted [15] which tries to improve on [13], with the introduction of a fully 
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3-momentum dependent mass. Here, the polarization is treated properly, but 
this leads to substantially more computational complexity, for no analytic 
form is found to approximate it in [15]. In these calculations r is found to 
increase with N, the number of flavours of fermion, with the value ranging 
from r 5.4 (for the N — 1 case) to r f=:i 14 (for the N — 3 case). The 
value given for iV = 1 is in good agreement with [14] , but the results of [14] 
suggest a slight reduction in r with increasing A^, not an increase. Inspection 
of the data of [15] suggests, however, that some caution may be necessary in 
accepting the results of the N = 2 and 3 cases. 

It is our view that the best way to proceed is by using the Matsubara 
formalism, due to the fact that the non-analj^icity problem does not arise. A 
simple closed form for the photon propagator to leading order in 1/A^ can be 
found [14]; and this leads to simpler equations than those found in the real 
time formalism, significantly reducing computational complexity. However, 
if one is to use the Matsubara formalism one must get round the problem 
of the i-r divergence. In [16] the divergence was regulated by introducing an 
arbitrary cutoff to take account of terms beyond leading order in 1/A^ in the 
vertex, which might regulate the equation for S. In this paper, we indeed see 
that by going beyond leading order in 1/A^ in the vertex we are able to get 
a finite equation forA^ — T,/A (with the inclusion of the transverse part of 
the photon propagator), where A{pf) is the wave function renormalization, 
without the introduction of an arbitrary cutoff. 

We shall divide the work in to four sections and an appendix. In the next 
section our interest lies in the formulation of equations for both '^{pf) and 
A{pf) at finite temperature for a completely general covariant gauge. If the 
gauge fixing function is a function of momentum, as it will be in Section.4, 
then this corresponds to a non-local gauge as discussed in [8] and [9]. The 
derivation of these equations depends on a simple choice of vertex ansatz 
for which T^{pf,kf) = A{kf)Y, which, although it does not satisfy the full 
Ward-Tahashi identities, is thought to be good starting approximation for 

In Section. 3, we go on to look at the Landau gauge. Here, we clearly see 
the problems associated with such a choice of gauge in the equations for S(pj) 
and A{pf); the poles in the kernals of both equations give rise to logarithmic 
divergences. Although this is the case, our next step is to see whether we 
can formulate a finite equation solely in terms of A4. Indeed, we find such 
an equation if we make the constant mass approximation Ai = Ai{0,7rT). 
However, if we choose our constant mass approximation to have the value 
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of Mipf) at some generally specified point p'^ — (p',pQy.), i.e M. — M.{p'f)^ 
when p' 7^ we see that the equation for M. contains divergences. 

In Section. 4 we show that for a choice of non-local gauge, which we shall 
call the D-gauge, the equations for A{pf ) and are both finite. Again 

we shall be able to derive an equation solely in terms oi M.{pf). By making 
the constant mass approximation Ai — M.{Q,t:T) we are able to show that 
this equation reduces to the same equation as in the Landau Gauge (making 
the same constant mass approximation), if we use the form of the photon 
propagator developed in [14]. However, there are two advantages to using 
the D-gauge. The first is that, unlike the Landau gauge, if we choose our 
constant mass to be = M.{p'f), where p'f is our generally specified point, 
we do not encounter any divergences. The upshot of this is that our equation 
for M. = Ai{0,7iT) can be considered a far more sensible approximation in 
this gauge. Furthermore, if one introduces full 3-momentum dependence into 
M{pf), the equation for M. remains finite. The second advantage is that at 
T — this gauge is identical to the Landau gauge. 

In the last section, Section. 5, we solve the equation for the constant mass 
approximation Ai = Ai{0,7rT). From our numerical solution we are able 
to calculate a value of r, which we find roughly to be r ~ 5.5. When we 
compare these results with those of [14], we see that these results compare 
favourably to those of [14] (using the constant mass approximation) with the 
number of flavours halved, the values of 2(7" = 0) being in good agreement. 

Finally in the appendix we prove an important result in our analysis. 

2. The Schwinger-Dyson equation 

For massless QEDs with A^-flavours we have the usual Lagrangian den- 
sity: 

c = -1/4/^.^'^ + J2Mi^ - (1) 

i 

where is the vector potential and i = 1,2...N. As in previous work ([14] 
and [16] ) we have chosen a reducible representation for the Dirac algebra, 
so that (1) has continuous chiral symmetry as discussed in [2]. The full 
Schwinger-Dyson (S-D) equation for the fermion propagator takes the form 
(in the Matsubara formalism) 

S-APf)=Sf-\pf)-- E y^5p(fc,)7^A,.(A;/-p/)r^(%-p/,M- 

(2) 
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In equation (2) we use (as in [14] and [16]) the subscript / to denote fermionic 
3-momenta, with kf = (A;o/,k); the zeroth component has the form fco/ = 
in the Matsubara formahsm, where m is an integer. To denote 
bosonic momenta we shall use the subscript b, with ph = {pob, p); the zeroth 
component is of the form pofo = where m is again integer. 

In this paper our interest lies in going beyond leading order in 1/N va. the 
vertex part. Our main motivation for this stems from the work of [14] and 
[16], in which it was shown that if the bare vertex is used, the transverse part 
of the photon propagator leads to an i-r divergence in the S-D equation at 
finite temperature. Furthermore, it was shown in [16], that the introduction 
of a dynamically generated fermion mass, E, into the fermion propagators in 
the calculation of the vacuum polarization, 11^^ did not alter this conclusion. 
The main hope of [16] was that if one went beyond leading order in 1/A^, 
such i-r divergences would be regulated, without the inclusion of an arbitrary 
cutoff (c.f [16]). 

Let us first write down a form for the fermion propagator, which takes 
account of wavefunction renormalization 

S-\kj) = A{kj)i> + nkf) (3) 

where A{kf) is the wavefunction renormalization and S is the fermion mass 
function. Now since we interested in going beyond leading order in 1/N in 
the vertex, we no longer take V to be its bare value of 67^^; instead we choose 
r*^ to be an ansatz of the form: 

T^{kf-pf.kf)=A{kf)eY. (4) 

Here one should note that although (3) may be considered as an exact ex- 
pression for S^^{kf), the expression for does not obey the Ward-Tahashi 
identities relating S~^{kf) to T'^{kf — pf,kf). Therefore (4) can only be 
considered at best an approximation of the true vertex. However, at zero 
temperature the vertex function is thought to contain only terms of either 
zeroth or first order in A(pf) (c.f [3] and [7]); our ansatz is thought to be a 
reasonable starting approximation to make for the vertex function. 

In our study we shall neglect terms beyond leading order in 1/A^ that 
contribute to the full photon propagator (these terms are thought to con- 
tribute little to the result (c.f [3]) at zero temperature); therefore we shall be 
using the approximate form for the full photon propagator (to leading order 
in 1/N) in the numerical calculations of the fermion mass, first developed in 
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[14]. For the moment, however, we shall write down a more general form for 
the full photon propagator: 



(5) 



where A^j, and B^y are the longitudinal and transverse projection operators, 
respectively. e(gb) is a general covariant gauge fixing function, which in gen- 
eral, we make dependent on 3-momentum. If e{qh) is not a constant w.r.t qi,, 
then one has a non-local gauge, as discussed in [8] at zero temperature. 

From (3), (4) and (5) we are able to write the S-D equation (fig.l) in the 
form 

S^iPf) ^ A{pf)i> + E(p;)/ = i^- J:{pf) (6) 

where 



g2 p (i^/c 



(3^ J {27rr ^fA{kf) + B{kf) [ 



DMA,. + DMB,. + !(|) 

(7) 



If we take the trace of (7), we get the following equation for 



Dr+DT + 



J 



(8) 



where = 

By multiplying by j) and taking the trace we are able to derive the fol- 
lowing expression for A{j)j): 



A{vs) = 1 + 



-2ri(po/, /Co/, P, k)D4gft) - 2r2(p, k)DT(g6) - 2^(p/.gft)(A;/.gft) 



where 71 = 



and Ti 



( (»-P/)8oA «i /, 



{(lb-kf)qob 



p^jk)B,, = k.p - 



k.qk.q 



(10) 
(11) 



Now that we have two equations (" gap equations ") for both E and A, 

the step is to investigate whether or not these equations are i-r finite, and 
to obtain an equation solely in terms of the physical mass, A4{pf), in the 
Landau gauge. This we investigate in the next section. 

3. The gap equations in the Landau gauge 

We shall first consider the Landau gauge, setting €{qb) — 0. The contribu- 
tion to the R.H.S of (8) most likely to contain i-r divergences is the po/ = ^o/ 
mode: 

^ <P^' = J J (2^ k^+pg, + A<^(k.po,) '^'-''i' + <^2' 

where -D°(q) and -Dj.(q) are the qob — contributions to Dl and Dt, re- 
spectively. Already, we know from the work of [14] and [16], that in the 

approximation A{pf) = 1, E = E(0,7rT) (which we call the constant mass 
approximation ), E° contains an i-r divergence at k^f = ttT. We now see 
whether a divergence exists for all pf at pf = kf. From the results of [10] 
and [14], we see that for small q 



where Mp = y^^f^- Already we see that any potential divergence must 
come from the contribution. By replacing by its small q behaviour 
and neglecting in (12), we have on performing the angular integration: 

_Ce^ r \k\d\k\ MjKpoj) 1 

P J (27r) k^+plf + M^{k,pof)\k^-p^\ ^ ^ 

where C is a constant of proportionality. One can see that (15) is manifestly 
divergent at k = p. An important point to note, is the appearance of a factor 
1//3 in front of the integral. This factor ensures that at T = there are no 
divergences in the equation for E, which is known to be finite at T = 0. 

Prom (15) and by comparison of (8) and (9) we see that A(pf) contains 
the same type of divergences. Although A{pf) and E(pj) both contain diver- 
gences, the question to ask is: can we find a finite equation for the physical 
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mass A4{pf), in which the divergences cancel? The answer is. indeed, that 
we can find an equation solely in terms of Ai{pf); however this equation 
generally contains divergences. 

To find an equation for 7W(p/), let us rewrite the L.H.S of (8) as A{pf )J^{pf) 
and then use equation (10) to eliminate A{pf), leaving us with 



P^J (27r) 



M{pf)pf{{pf.kf){DL{qb) + DT{qb)) - 2T^{Pof, P, k)Dj.(g,) - 2T2(p, k)DT(gb)) 



M{kf){DL{qb) + DT{qb)) 

.(16) 



k) + M\kf) + k) 



k} + M^{kf) + k} 

Now, our first goal would be to solve this equation for a constant mass 
A4 = M.{pf). Instead of specifying the physical mass at the point (0, ttT) 
for all pf (as is usually done in similar calculations), let us choose A4 to be 
the same as the value of A1(p/), specified at some general point = (pg^, p') 
for all Pf. In this approximation, M. — M.{p'f), equation (16) becomes 



Dr^ikf - p'^) + Drikf - p'f) - l^[{p'^.kf){D,^{kf - p'^) 



(3 ^ J {2tt)^ + k} [ ' ' pf 

+DT{kf - p'f)) - 2TMf, kof, p', k)DLikf - p'f) - 2T2(p', \^)DT{kf - p'f)]] .(17) 

Possible danger comes from the kf = p'j pole in (17), so we look at the 
contribution to the R.H.S of (17) where kof = pg^, which we call A^°. 

(p-..- ^--"--/2;jr'''^ p°)^^(''-po). (18) 

Terms in M.^ proportional to -D°(k — p') are finite, for they are regulated 
by Afp in the denominator (c.f. (14)). We are also able to show that the 
contribution from 

is finite (see Appendix) , although at k = p' there exists a term proportional 
to a delta function which vanishes as p' — > 0. This leaves us with a term 
proportional to p'^L'j,(k — p), which is divergent when we choose p' ^ 0. 
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The upshot of the above is that although we have been able to find a 
constant mass equation which is finite at p' = 0, this approximation is not a 
sensible one to make, due to the wild discrepancies between points where p' = 
and those where p' 7^ (the equation for p' 7^ leads to a singular mass). 
These singularities are also manifest if one decides to relax the constant mass 
approximation. In the next section we shall work in a gauge which, as we 
shall see, removes these singularities. 

4. The D-gauge 

If we choose a non-local gauge £(^5) = {D^ — DT)ql (which we call the 
D-gauge) then, already, we see from (8) that 



S(P/) 



2e2 



E 



(fk Mjkf) 
{27ryk} + M(kf) 



DM, 



(19) 



which is finite for all pj. In this gauge we are also able to rewrite (9) as 



1 



E 



pjP^J {27rr kj + M^ikf) 



[2{pf.kf)DL{q,) 



-2Ti(po/, kof, p, k)Dr.iq,) - 2T,ip, k)L>^(5,) + ^bz-g^Kfe/.g^) ^^^^^^^ _ ^^^^^^^ 



.(20) 



Again the contribution most likely to contain divergences in the R.H.S of 
(20) is the kgf — Pof contribution. 



x2 



(p.k) 



(p.q)(k.q 



1 



d^k 



r d^l 
J (2^ 



pjpj i27:rk^+plf + M{bfk,pof) 
2(p.q)(k.q 



{p.k) 



D° (q) . (21) 



Since we already know that the term 

2k.(k-p')p.(k-p') 



p.k 



(k - p)^ 



D^k - p') 



gives a finite contribution (see appendix), it is easy to see that (21) is i-r finite. 
The consequence of this is that we have a gauge in which our equations for 
E and A are finite. An advantage of this gauge (as opposed to other gauges 
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we might choose to cancel this divergence) is that in the hmit T ^ 0, this 
gauge joins smoothly on to the Landau gauge, for at T = = Dt- So at 
T = 0, this gauge has exactly the same gap equations as those of the Landau 
gauge, which has been extensively studied. 

Prom (19) and (20) we are able to derive an equation solely in terms of 
the physical mass, M.: 



2M{kf)DL{qb) 



Mipf) 



[2{pf.kf)DL{q,) - 2Ti(po/, fco/, P, ^)DL{qb) 



-2T,(p,k)D,(,.) + ^(^/•^''f/-^^) (!,,(,,) _ D.iq,)) 



(22) 



This equation, unlike (16), is finite for all p/, due to the fact that both (19) 
and (20) are finite. Again, as in the previous section, we make the constant 
mass approximation, now choosing = (0,7rT). In this approximation (22) 
reduces to 

-2 (Pk 



-2{'KTf 1 - 



k 



06 



kl 



DL{h)\ + 



(27r) 
27iTkQbkf.kh 



2DL{h) 
M^ + kj 



;{2nTkof)DL{h) 



'-{Drikb) - DL{kb)) 



1 



M^ + k}^ 



Now using the approximation for the photon propagator given in [14], (23) 
becomes: 

d^k 



4D^{kb) 2kof D^ih) 



(2vr) 



2k^ob D^{h 



^M^ + kj ttT M'^ + k} 
(Pk 2Dl{h) 



a 



kl M^ + k}^ 



N(3J {2tiY M^ + k) 



(24) 



where a = Ne^ 



kl + ah/9> 
1 

k^Tnp) 



^06 7^ 0, 



k 



06 



0, 



andnO(k)4g^ 



,^ 161n2 (-Ti(3\k\ 
|k|p H — exp 
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(25) 
(26) 
(27) 



One can easily check that, in the Landau gauge, using the same constant 
mass approximation ( Ai = A^(j9qj,p') = 7Vl(7rT, 0)) and the same form 
for the full photon propagator, (17) leads to the same equation as (24). 
However, in this non-local gauge (24) seems a more sensible equation to 
solve, due to (22) having no divergences for any pf. So if we wanted to solve 
(22) with a constant physical mass approximation, specified at some other 
point, say A4 = Ai{p'j-), we would still get a finite result for A^. This gauge 
has an added advantage; if we wanted to solve (22) with full 3-momentum 
dependence, we would still get sensible results. 

5. The numerical solution to the constant physical mass equation 

(24) 

By following similar steps to those of [14], we are able to recast (24) in a 
form convenient for numerical analysis 

1 = — V [(1 - 2n)/i(27rn, 0.125a, (aV + (27r(n + 1/2))^^^) 



n=l 



+ (1 + 2n)/i(27rn, 0.125a, (aV + {2n{n - 1/2))^)^/^ 
-(27rn)2(/2(27rn, 0.125a, (aV + (27r(n + l/2)fY/^ 
+I2{2nn, 0.125a, (a^s^ + (27r(n + l/2) f y^^) 

X 1 



+ 



a r 
dVio 



dx- 



ttN Jo '^x^ + f3^Ul{x)x^ + a'^s^ + n'^ 
where s = M./a, a = a(3 and x = /3|k|; also 



(28) 



Ii{d, a, c) 
and 

l2{d,a,c) 



1 



In 



a 



2(a2 + c2 - ^2) y(^d + a)^J (c2 - ^2)1/2(^2 + ^2 - ^2) 



arctan 



(29) 



a 



d^)d 



+ 



(a^ + - d'^)\nd - {c^ - d^) \n{d + a) - aMn c' 
a2(c2 -d2)(a2 + c2 -(i2) ^ 

a f(c''-d^Y/^\ 



(c2-d2)i/2^ 



(a2 + c2-d2)(c2-d2)3/2 



arctan 



.(30) 



Before proceeding to solve (28) numerically we require the zero temperature 
limit of (24). This is most easily found by going back to (22), due to the 
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fact that the sums in (28) and (24) cannot be performed analytically. In the 
T ^ limit (22) becomes: 



2p'^M{k) 2{p.q){q.k) M{p) 



M^{k) + A;2 



(31) 

By Taylor expanding (31) in powers of 3-momentum, p, keeping only terms 
up to order and performing the angular integration, we are able to get the 
following expression for s — A4/a in the p — > limit: 



.(0) 



(27r)2iV 



dy- 



TT: / dy- 



s{0) 



lo {y + 0.125) {y^ + s^{y)) 12 Jo {y + 0.125)^ {y^ + s^{y)) 

(32) 

where y — k/a. Since S{y) is monotonically decreasing, so as to satisfy (22), 
it is true in general that 



s{0)> 



{27v)m 



/ dy 
Jo 



1 - 



(y + 0.125){y^ + s^{y)) 
\ -1 



dy 



1 



6(27r)2iVio (y + 0.125)2 + (y2 + s(o))^ 



(33) 



If we look at (33) we see that the effect of wavefunction renormalization is to 
increase s(0), for if we set A = 1 we neglect the second term in the R.H.S of 
(32). In the constant mass approximation s{y) = s(0) this bound becomes 
saturated: (33) becomes an equation for s(0). On performing the integrals 
in the constant mass approximation (33) reduces to: 



where Fi 



(0.125)^ 



TT 



(27r)2Ar 
0.125 



and Fr 



(0.125)2 + ^2^25 (0.125)2 + s2 
/ (0.125)0.25 0.25(0.125)2 



[F^{s) + F,{s)/3] 
In 



(0.125)2 + ^2 ((0.125)2 + s2)2 



0.125 

TT 

2^ 



+ 



0.125 



(0.125)20.25 
(0.125)2 + ^2 ~ ((0.125)2 + s2)^ 



In 



.0.125 
0.125 



(0.125)2 + s2 



(34) 
(35) 



(36) 
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Now we are able to plot graphs of the solutions as functions of (kBT/a) for 
various (see fig. 2), by solving (28) and (34) numerically using Mathemat- 
ica. If one halves the number of fiavours given for each solution in fig.l, 
these results differ little from those of [14], the values of E being in good 
agreement. In fig. 3 we show a graph of TJzbIol^ the critical temperature 
plotted against N. As expected, we see that as N increases TJzb falls. How- 
ever, there is no finite A*",,, a critical number of fiavours. This is due to the 
singular behaviour of the R.H.S of (34) at s = 0, from which we infer that 
s = is not a solution of (34) at a finite value of N\ therefore iVc — >■ oo, 
as was shown in [14] and [16] for the A — — T^{7tT, 0) approximation. 
There is strong evidence to support the existence of a finite Nc in the full, 
untruncated mass gap equations in the absence of any approximations. The 
reason for this discrepancy is that our approximations are cruder than those 
used at T = 0. Even though this is the case, we expect the exact solution 
(at the point (ttT, 0)) to have, at least, the same qualitative features as those 
shown in fig. 2. 

Finally we present a table of values (Table. 1) for r = which changes 
little with and takes the average value r = 5.56. When compared to the 
results of [13] and [14], it is evident that wavefunction renormalization does 
httle to change the value of r. The same conclusions were found when [12] 
was compared with [11] in the unretarded case. 

An obvious extension to this work would be to relax the constant mass 
approximation by introducing frequency or momentum dependence (or both) 
into our equation for the physical mass. Another possible extension would 
be to use a more sophisticated vertex ansatz. A good one to choose would 
be the Ball-Chu ansatz which obeys the Ward identities at zero temperature. 
The main difficulty is that the form of the equations for A and S would be 
far more complicated than those encountered in this paper. As this may 
well be, one long term goal should be to use an ansatz which respects the 
Ward identities at finite temperature, which would reduce to the Ball-Chu 
ansatz at zero temperature (the Ball-Chu ansatz does not respect the Ward 
identities at finite temperature, due to loss of Lorenz invariance). Such an 
ansatz would get rid of the divergences seen in the Landau gauge and would 
produce a gauge dependent physical mass. Until then the (L'l-Dt)?^ 

gauge 

is a sensible choice of gauge to work in for two reasons already stated: firstly 
this equation leads to finite equations for A{pf) and S(p/), secondly at T = 
this gauge is identical to the Landau gauge. 

From this paper, one can plainly see that there is still a considerable 
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amount of work to be done at finite temperature in QED^, especially if we 
want our results to be consistent with new results at zero temperature ([7], [8] 
and [9]). 

Appendix 

In the appendix we shall show that the term 

2k.qq.p 



p.k 



^^(q^k-p) 



(A.l) 



in (21) leads to a finite integral when k — p — > 0, although we shall see that 
a delta function contribution exists at k — p = 0. To show this we consider 
the integral 



27r 



q2 / q2 



In (A.2) we have replaced D^(k — p) by its leading order behaviour in q. 
This is sufficient, for wc are interested whether or not there is any divergent 
behaviour at p = k; higher order corrections do not matter, for they are 
unable to cause divergent behaviour when integrated over k (the next highest 
term is proportional to l/|q|). We are able to write (A.2) in terms of a 
complex integral (where we have introduced a regulator e, so that we are 
able to probe the k = p region.) 

1 



27ri J\z\=i 



2p2k2z- (k2 + p2)|k| 


IpI 


(.2 + i)/2 




k2p2(-S- 


p|/|k|+Z6/(2|k||p|))2(^ 




k|/|p 


-Z6/(2|k||p|))2 



(A.3) 



We are able to evaluate this integral, giving us the following expressions. 



7 = 



2p'k' - {k' + p2)(p2 - ie/2) (4|pHk| - 22|p||k|e)|p||k| 



(p2 — k^ — ie)^ 
(k^ + p^)kV(pVk2 - ze/k^ - eV(4kV) + 1) 
(p2 — k^ — ie)^ 



(p2 - k2 - ie)3 
when IpI < |k|(l-52)VfA.4) 



7 = when |k|(l + S^^^ > |p| > |k|(l + S^f^-^) 
2p2k2 - (k^ + p2)(k2 + ie/2) (4|k|3|p| + 2i|p||k|e)|p||k| 



(k2 - p2 + ie)2 

(k2 + p2)k2p2(p2/k2 - te/k" - e2/(4k2p2) + l) 

(p2 — k2 — ie)3 



(k2 - p2 + ie)3 
when IpI > |k|(l + (52)^fA.6) 



13 



where 5 — 2\k\\p] - checking (A. 4) and (A. 6) it is easy to show that the 

terms of zero order in e, in the numerator, vanish in both expressions when 
e — s> 0. For terms of first order in e (in the numerator) we have the following 
expressions 



-ie/2 



eV2 



(p2 — — ie) 



+ie/2 



+ 



- k^ — ie)"^ (p^ — k^ — ie) 

when IpI < |k|(l- (5^)1/2 
eV2 kh^ 



(k2 - p2 + ie) (k2 - p2 + ie)2 



+ 



p'^ — k^ + ie) 
when IpI > |k|(l + 52)V2_ 



(A.7) 



(A.8) 



In both (A. 6) and (A.7) the first two terms cancel when e ^ 0, but the last 
terms in both expressions is singular at k^ — p^ = (zero everywhere else). 
When we consider terms that are second order in e (in the numerator) we 
see, that also, these terms are singular at k^ — p^ = 0. 



iil±i?4^when|p|<|k|(l-^^)V^ 



(^p2 — k^ — ie)^ 

,2 I „2\^2 



7,2 



(k2 + p2)e74 , , , „ , ,^ .2M/2 

-^^ / when p < k 1 - 5'^/^ 

(k^ — p"^ + te)"^ 



(A.9) 
(A.IO) 



Combining (A.7) with (A.9) and (A.8) with (A.IO) gives us for the singular 
part of /: 



"■sing 



[p2 _ + p2)/4] 



'■smg 



/„„5=[k'-(k' + pV41 



[j^TZ^r-^) when|p|<|k|(l-^Y'^ 
when |k|(l + 8^f/^ > IpI > |k|(l - 5^f/^ 

when IpI > |k|(l + 0^/fA.ll) 



(k2 -p2 + ie)2 



In the hmit e — > 0, (A.IO) suggests that Igmg will be proportional to a delta 
function 



' smg 



[V^/2)M5{\^ - p2) 



(A.12) 



where we are able to evaluate N 



H = hm 



e^du 



f^o [-/s^k^ [u + ie)^ 
14 



eHu 



[u + ie)' 



6^-^ \ 2((52k2 + ie)2 + 2{ie - S^k^)^ J 
= -1. (A.13) 
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r 


1 


5.59 


2 


5.57 


3 


5.51 



Table 1: The value of the quanity r as defined in the text at the indicated 
values of N. 

Table 1 

Figure Captions 

Figure. 1: A diagramic representation of the S.D equation. 

Figure. 2: Graph showing the constant physical mass solutions as func- 
tions of T/cb/q! for A'" = 1, 2 and 3. 

Figure.3: Graph showing the variation of ksTc/a with w.r.t N . 
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Fig. 1 
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Fig. 3 
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